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Transmission of an arbitrarily poiarizid plane wave by an arbitrarily oriented spheroid in the 
short- wavelength limit is considered b y-tH e -g a e of ray theory. The transmitted electric fieid is added to 
the diffracted plus reflected ray-theory electric field that was previously derived to obtain an approxima-</^ fl 
tion to. the far-zone scattered intensity in the forward hemisphere. Two diiTerent types of cross- 
polarization effects are found. These are (a) a rotation of the polarization state of the transmitted rays 
from when they are referenced with respect to their entrance into the spheroid to when they are 
referenced with respect to their exit from it and (b) a rotation of the polarization state of the transmitted 
rays when they are referenced with respect to the polarization state of the diffracted plus reflected 
rays. 0 1996 Optical Society of America 


1. Introduction 
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A number of optical particle-sizing instruments, such 
as the phase-Doppler particle analyzer, measure the 
diameter of small particles by analyzing certain fea- 
tures of their forward-hemisphere light-scattering 
signature. 1 * 2 For particles much larger than the 
wavelength of light, forward-hemisphere scattering is 
well approximated by ray theory, which considers the 
diffracted, reflected, and transmitted rays that reach 
the detector. For scattering by a sphere, the ray- 

theory far-zone scattered intensity closely matches 

A'X > the results of Lorenz-Mie theory 3 * 4 for particle-size 
parameters as low as 30 and for scattering angles as 
large as 50°. As a result of this close match, phase- 
Doppler particle analyzer calibration curves have 
been based almost entirely on ray theory. 1 * 2 * 5 * 5 
A basic assumption underlying the calculation of 
optical sizing-instrument calibration curves is that 
the particles being sized are spherical. A necessary 
first step in evaluating the response of these instru- 
ments to nonsphencal particles is to nave an accurate 
and easily implemented theory of light scattering by 
such particles. Because the simplest nonsphencal 
particle is a spheroid, understanding scattering by a 
spheroid should be helpful in determining the impor- 
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tant features of scattering by nonsphencal particles 
in general. 

There are a number of different methods for ex- 
actly solving the electromagnetic boundary- value prob- 
lem of a plane wave scattered by a spheroidal par- 
ticle. 7-11 The numerical implementation of these 
methods, however, suffers from ill conditioning for 
size parameters greater than approximately 35 and 
for large spheroid eccentricities. 3 * 11 * 12 The rav- 
theory model of spheroid scattering is an attractive 
alternative because it is expected to provide a reason-s 
able approximation to the solution of the exact wave- 
scattering problem for particle-size parameters begin-' 
ning approximately where the numerical imple-- 
mentation of the exact methods starts becoming ill 
conditioned. Ray scattering by a spheroid has al- 
ready been applied to the analysis of the generalized 
rainbow caustic in the backward hemisphere, which 
is caused by the confluence of a number of rays 
making one internal reflection within the spheroid 
beiore exiting. 13 ” 19 Ray-tracingprograms for scatter- 
ing by an arbitrarily shaped particle have also been 
developed for certain specialized applications. 20 - 21 
The purpose of this paper, along with a companion 
paper 22 that is hereafter designated as part I, is to 
describe scattering in the forward hemisphere of an 
arbitrarily polarized plane wave by an arbitrarily 
oriented dielectric spheroid by the use of ray theory. 
In part I diffraction and specular reflection were 
considered. In this paper transmission, which is the 
third of the physical processes expected to dominate 
scattering in the forward hemisphere, is considered. 
Diffraction plus reflection was exactly soluble in the 
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sense that the magnitude, phase, and polarization of 
the electric field of the outgoing rays were directly 
> expressible in terms of the scattering angles 0 and <t>. 
This is not the case for transmission because of the 
complexity of both the refraction geometry and the 
spheroid shape. If the plane wave is not incident 
parallel to the spheroid major axis, another complex- 
ity occurs as well. For this case the plane of inci- 
dence of a ray at its point of entrance on the lit side of 
the spheroid does not coincide with the plane of 
incidence at its point of exit. This leads to cross- 
polarization effects that do not occur for scattering by 
a dielectric sphere. 23-25 For example, a transmitted 
ray that was incident upon the spheroid with the 
transverse electric (TE) polarization will exit it ellipti- 
cally polarized, i.e., with a mixture of TE and trans- 
verse magnetic (TM) polarizations. Similarly, the 
TE polarization directions of the transmitted and 
reflected rays reaching an observer at the scattering 
angles 0, <t> are rotated* with respect to each other. 
The transmitted ray is additionally eiliptically polar- 
ized when referenced with respect to the TE and TM 
polarization directions of the reflected ray. 

Cross-polarization effects also occur for scattering 
by a nonspherical particle much smaller than the 
wavelength of light because of differences in the 
particle's polarizability in different directions. 26 - 27 
Similar effects also occur for an optically active 
particle. 28 Cross-polarization effects originating 
solely from the geometry of the particle also occur for 
scattering by a dielectric cylinder at diagonal inci- 
dence. 29 For scattering by a spheroid, the cross- 
polarization intensity has been computed with the 
exact solution to the wave-scattering problem. 30 
Our purpose here is to demonstrate clearly and 
explicitly the geometric origin of the cross-polariza- 
tion effects by-feke-trse-itt the context of ray theory. 
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ne b(Jdy“5T this paper proceeds as follows. Sec- 
tion 2 presents a brief review of both the notation and 
the spheroid geometry that were described in detail in 
part I. In Subsections 3.A.-3.C. the magnitude, 
phase, and polarization vector of the transmitted 
electric field are obtained. In Section 4 these results, 
are combined with the diffracted plus reflected elec- 
tric field obtained in part I, and the cross-polarization 
effects are examined. Last, in Section 5 the n uiwer i - 
sal computation of the diffracted plus reflected plus 
transmitted far-zone intensity and compa re-our re- 
sults with those obtained by the\ise of other methods. 
( / 
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2. Spheroid Geometry 

Consider a plane wave of wavelength wave number 


9— 


(i) 


and angular frequency u> propagating along the posi- 
tive z axis of an xyz coordinate system fixed in the lab 
reference frame. The electric field of the plane wave 
is 


£K 

E inc = £ 0 (cos x“z + sin x« y )exp (ikz - iwf). (2) 


The plane wave is scattered by an arbitrarily oriented 
spheroid. The transmitted contribution to the out- 
going electric field in the far zone at the scattering 
angles 0, <P is of the form 


Cl< 

<Ail 


EtraJ©. <*>) = ££ exp (ikR - iuit)S lnta (Q, <t>) 

x exp[i5 trsns (6, tytaaje, <t>). (3) 

C K 

In Eq. (3) R is the distance from the center of the 
spheroid at the origin of coordinates to the observer; 

S trans is the magnitude of the transmitted electric 
field, ignoring the Fresnel coefficients at the spheroid 
surface; 5 trans is the phase of the transmitted electric 
field with respect to that of the reference ray of 
Subsection 3.D. of part I; and is the polarization 
vector of the transmitted electric field and contains 
the Fresnel coefficients. The dielectric spheroid has 
the real refractive index n > 1. There is no other 
restriction on the numerical value of a. The surface 
of the spheroid is given by 



(4) 


where the xy'z axes are attached to the spheroid. 
The spheroid z" axis lies in the 9, 4) direction with 
respect to the xyz lab coordinate system, where 0 ^ 

9 ^ tt/ 2 and 0 ^ <J> < 2 t t. The description of. 
scatteringby the spheroid is simplified considerably if 
we use a second lab coordinate system x'yz rotated 
with respect to the .ryz lab system by the angle 4> 
about the z - z f axis. In this new coordinate system 
the equation of the lit (i.e., lower) and the shadowed 
(i.e., upper) spheroid surface is 


T 

-upper 

-lower 


= uuAr cos 



/2U/2 


(5) 


where 


lb 2 - a 2 \ 


w - sin 9 cos 91 1 • 

(6) 

A = ( b 2 sin 2 9 + a 2 cos 2 0) I/2 , 

(7) 

B - a, 

(8) 


and where the elliptical coordinate r and are 
defined by 

x' = Ar cos 5', 

y' = Br sin (9) 

with 0 £ r ^ 1 and 0 ^ < 2-r. The outward unit 

normal to a point on the lit surface of the spheroid, 
m , and the outward unit normal to a point on the 
shadowed surface, n , are shown in figure 4 of part I, 
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Fig. 1. Geometry of the transmitted ray. The unit normal m 0 ' to 
the surface at the point of entrance is in the t - y 0# -n 0 direction 
► with respect to the x\ y' , z' rotated iao coordinate system. The 
portion of the ray inside the spheroid has length s Ql and is in the 
^ 01 . Tloi direction. The unit normal /if to the surface at the point 
of exit on the shadowed side of the spheroid is in the r| t 
direction. 


and their equations are given by equations 14 and 15 
of part I. 


Using Eq. (14) ofpart I and Eq. (5) above, we obtain 


0 ^ 


ab 

“yj Qo sin § 0 ' 

tan no = ^ . (12) 

a o' cos <j 0 ' + w 

A~ 


tan 3 (— - ^ 0 ) = 


jab 

1^5 <7o' cos ? 0 ' 



with 


f ab t \2 
^ p sin I , 

(13) 


<?o = 


ro 


(1 - r 0 ' 2 ) 1/2 


(14) 


These results are analogous to equations 21, 23, and 
24 of part I for the normal to the lit surface at the 
point of incidence of the reflected ray. 

After transmission into the spheroid, the ray propa- 
gates in the 'f'm, not direction with respect to the 
x'y'z axes. The unit transmitted wave vector is 
given by 


£ t0 = sin '('a! cos "Hoid*' 
+ cos ^oid r '. 

It may also be written as 31 


sin sin T| 01 d 


(15) 



cos 0 iO 


— cos 0 fO hn 0 ' 


where 


(16) 



A. Magnitude of the T ransmitted Electric Field 

The trajectory of a ray transmitted through the 
arbitrarily oriented spheroid is shown in Fig. 1. All 
quantities pertaining to the ray as it enters the 
spheroid on the lit surface have the subscript 0, all 
quantities pertaining to the ray inside the spheroid 
have the subscripts 01, and all quantities pertaining 
to the ray’s exit from the spheroid have the subscript 
1. An incident ray is parameterized by the rotated 
lab frame coordinates and (jo’. The unit wave 
vector of an incident ray is 


*.o = u,'. (10) 

The unit normal rh 0 ' to the lit surface at the point of 
entrance of the ray has the spherical angles 77 — 'I'q, no 
with respect to the.t'y'z' axes, i.e., 


"V- cos n 0 “ r ' 4- sin '(' 0 sin n 0 «/ - cos 

( 11 ) 


0 ;o = 'I' 0 (17) 

is the angle of incidence of the ray at the lit surface. 
The angle of refraction 0, o is given by Snell’s law 

n sin 0 (O = sin 4' 0 . (18) 

A comparison of Eqs. (15) and (16) allows us to relate 
the angles of h t Q to the angles of itiq ' . We obtain 

*01 = ~ 9(0 

hd = fO + - (19) 

The relation between not and no is expected on 
physical grounds. The plane ot incidence containing 
k i0 and rn.Q makes an angle no with respect to the x 
axis, so n,o must also lie in this plane. Furthermore, 
Fig. 1 suggests that dig and k,o lie on opposite sides of 
the z axis in the no plane accounting for the differ- 
ence of — between no and not- 
The trajectory of the ray inside the spheroid is 
given by 
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*' = V - Sol sin * 0l cos t]q, 

y' = yo' ~ s 01 sin * 01 sin t,o. 

= V + S 0l cos * 0l , (20) 

where Sg, > 0. The transmitted ray exits the spher- 
oid on either its shadowed or its lit surface. By 
combining Eqs. (5) and (20) we find that the distance 
of travel of the ray inside the spheroid for either 
possibility is 


s 0l (cos * 0l + w sin * 0l cos no) - — (1 - r 0 ' 2 ) l/2 

Thus the sign of the left-hand side of Eqs. (25) and 
(26) may be used as an indicator as to which side of 
the spheroid the transmitted ray exists. 


2ab 

s oi = — (1 ~ r 0 ') t/2 


. ab I B 

(cos * 0I -r w sin * 0l cos no) + Jg ?o' sin * 01 sin § 0 ' sin no + 7 cos cos no 


9i 9 

a-b~ 


(cos * 0 i + w sin *01 cos no) 2 + 77^5 sin 2 * 01 sin 2 t| 0 + — cos 2 


A 2 B 2 


E l 


A 2 


>T lo 


( 21 ) 


The coordinates of the ray at its point of exit from the 
spheroid are found by combining Eqs. (20) and (21) to 
give 


z i' ~ x o' ~ s oi sin * 01 cos no * Ar,' cos 
y i' = yo' ~ Sqi sin * 01 sin no = Br x ' sin gj', 


= *0 + s 0 i cos *oi. 
which is equivalent to 
„ 2r 0 ' 

r i ‘ = r o ' - -g- s 0l sm *ot 

B 


( 22 ) 


x sin £ 0 ' sin no + T cos & cos no 
\ A 


s oT . 9 „ B 2 . 

^ 2 ^ ^ m)i sin no *** "^9 cos- t| q | t 


( 23 ) 


r o sin § 0 ' - — sin * 0l sin no 
tan = • (24) 

, _ , S Q\ . 

r o cos $ 0 - “T* sin cos 


Let us assume for the moment that the ray exits on 
the upper or shadowed side. The unit normal n , ' to 
the spheroid at the point of exit has the spherical 
angles * tl m with respect to thejc'yV axes, i.e., 

ti l = sin *, cos m^x' + sin *, sin ni^y + cos *,i2.', 

(27) 

Using equation lo of part I and Eq. (5) above, we then 
obtain 

ab 

AB qi ' sin ^ 1 ' 

tan n i = ^ (28) 

^<7i' cos - w . 

tan- *i = q{ cos - w) + ?i’ sin §,'j , 

( 29 ) 

, r [ 

“ M - r '2U/2 * (30) 


For most commonly encountered values of 6/a and 
n, the transmitted ray exits the spheroid on its 
shadowed side. For this case, combining Eqs. (5) 
and (22) gives 

, r ab 

s o\( cos '[' 0l -r W sin 'I' 01 cos Tj a ) - — (1 - r 0 ' 2 )V- 

ab 

= jd -rr 2 ) ,/2 . (25) 

But for a spheroid with high eccentricity and large 
refractive index, the point ol exit of the transmitted 
ray sometimes occurs on the lit surface. For this 
case, combining Eqs. (5) and (22) gives 


We next determine the angles of incidence and 
refraction of the ray at its point of exit. The incidenc- 
unit wave vector at the point of exit is 

~ k a = k<o- ( 31 ) 

Its explicit form is given by Eqs. (15) and (19). The 
angle of incidence 9 a for exit on the upper surface is 
given by 

cos 0 :1 = h/ * k ti 

= cos * 0l cos *t - sin * 0l sin *, cos(-n 0 - q,). 

(32) 

The transmitted angle 9 ( , is again given by Snell’s law 
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sin 0 fL = n sin 0 {1 . (33) 

If n sin 9 a > 1, the ray is incident upon this interface 
past the critical angle for total internal reflection and 
transmission does not occur in ray theory. In actual- 
ity, light waves are transmitted and are described by a 
Fock transition. This effect is not modeled in this 
paper. 

For n sin 0 fl ^ 1, the ray is transmitted out of the 
spheroid and its final outgoing unit wave vector is 

k a — sin 0 cos( ( b — <f>)u x ' + sin 0 sin(<$ - <fr)u y 

+ cos 0z2 s ', ^4) 

which may also be written as 

k n = nkn + (cos 0 (l - n cos ' (35) 

Comparing Eqs. (34) and (35), we find that the 
scattering angles 0, are given by 

cos © 

= n cos 4^ 01 + (cos 0 fl - n cos Q^cos (36) 

tan <t> 

(cos 0 fl ~ n cos fyjsin jh sin t| 1 — n sin j'pi sinn 0 

= (cos 0 a - n cos 0 a )sin cos - n sin ^ 01 cos r\ Q 

(37) 

Thus given the coordinates and §o' of an incident 
ray, its scattering angles © and <D are obtained by 
combining Eqs. (12) — (14), (18) and (19), (21), (23) and 
(24), (28H30), (32) and (33), and (36) and (37). 

If the transmitted ray exits on the lower or lit side 
of the spheroid, the unit normal rhy is taken to have 
the spherical angles ~ “Hi vnth respect to the 

xy'z axes. For this case, the factors of -w in Eqs. 
(23) and (29) should be replaced by +u/, and the 
factors of cos Th in Eqs. (32) and (36) should be 
replaced by -cos '{'V 

These expressions for the scattering angles are 
complicated because of the cumbersomeness of the 
spheroid geometry and the cumbersomeness of the^ 
geometry of refraction, and because the orientation ot 
the spheroid is arbitrary. For arbitrary incidence it 
is easily shown that the planes of incidence at the 
ray’s point of entrance into and exit from the spher- 
oid do not coincide. Specifically, the first plane of 
incidence makes an angles n 0 with respect to the x* 
axis. But the second surface normal makes an 
angle r\ x with respect to the x' axis. If t) Q * r) lf the 
two planes of incidence are not coincident. 

If rj 0 = r[ U the two planes coincide and the entire 
trajectory of any given ray lies in a single plane. 
Comparing Eqs. (12), (24), and (2S), one can see that 
this occurs only for w = 0 and A — B corresponding to 
end-on incidence, where the z" spheroid major axis 
lies in the propagation direction of the incident plane 
wave. For the end-on geometry 

= flo = m' = = tf> - <i> _ "■ ( 3S ) 


Equation (13) then simplifies to 


tan^o = -<7o'. 


Eq. (21) simplifies to 


cos 'Koi + ” qo' sin 

Soi = 26(1 - r 0 ' 2 ) l/2 | a — 2 ’ (40) 

cos 2 'f'o, + — sin 2 
a~ 




Eq. (23) simplifies to 

f 2r 0 's 0 i . 


sin x V Ql + — sin 2 'V Qi 
a~ 


, sox . 

= r 0 - — sin V F oi , 


Eq. (29) simplifies to 


6 

tan 'P 1 = ~Qi> 

(42) 

Eq. (32) simplifies to 


8il = *01 + *1. 

(43) 

and Eq. (36) simplifies to 


> 

1 

CD 

II 

CD 

(44) 


Unfortunately, Eqs. (19) and (39)-{44) are still suffi- 
ciently complicated that they cannot be analytically 
inverted to obtain r 0 ' as a function of 0. 

The inability to obtain r 0 ' and & as functions of 0 
and <l> prevents us from analytically evaluating the 
magnitude of the transmitted electric field. 


S^(Q, *) = 


krABr Q ' \ 1/2 
d sin 0/ * 


in terms of 0 and <f> alone, where 
0 © 00 

4 . ^ W . (46) 

dr a ' d£, 0 ' CK_ 

(See equations 17 and 18 of part I). In Section 4 the < (\ 
numerical evaluation of Eqs. (45) and (46) is dis- 
cussed. 

B. Phase of the Transmitted Electric Field 

The trajectories of both the transmitted ray and the T~ 
reference ray are shown in Fig. 2. The reference ray if . — 
is described in detail in part I. Briefly, it propagates 
along the c axis to the origin as if the spheroid were 
absent. It then turns to the 0, T direction and 
propagates in that direction to the far zone. The 
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Fig. 2. Trajectories of the transmitted ray and the reference 
ray. The spheroid entrance plane is UU\ and the spheroid exit 
plane is W\ 


distance a/? from the spheroid entrance plane UU' to 
the origin is 


~min ► (4 0 

where z min ' is the lowest point on the spheroid 
surface. The distance from the origin to the spher- 
oid exit plane W' is (equation 41 of part I) 


3/? - 


A 2 


— [cos 0 - w sin 0 cos(<t> - <f ))] 2 " V A 


a 2 b 


sin 2 0 


cos 2 (<(> - 4>) sin 2 (<ft - 4?) 


A 2 


B 2 


- (48) 


The optical path length of the transmitted ray 
between the spheroid entrance and exit planes is oq *f 
fis 0 i + p t , where a 0 is the distance from the entrance 
plane to the ray’s point of entrance on the spheroid's 
lit surface and p t is the distance from the ray's point 
of exit on either the spheroid's shadowed or lit surface 
to the exit plane. From Fig. 2, the distance a 0 is 
given by 


«o = + V( r o\ So')- (49) 

The distance 0 t is determined as it was for reflection 
in part I. The equation of the outgoing reference ray 
is 

x' = (3 ff sin 0 cos(<f> - <(>), 

y* = p/? sin © sin(<b - 6), 

z' - 3/? cos 0 (50) 

for > 0. The equation of the plane normal to this 
line at the point 0* given by Eq. (48) may now be 
constructed. This is the exit plane W' of Fig. 2. 


The intersection of this plane with the outgoing 
transmitted ray 

x* = Ar { f cos + 0! sin 0 cos (<t> - <b), 

y* = Sr L ' sin -r 0 : sin 0 sin(<f> - <b), 

2 ' = 2i'('Y, Si') + .0! COS 9 (51) 

for 0t > 0 is then determined, giving 

Pi = Pa - Ar L ' cos sin 0 cos($ - 4>) 

- Br{ sin sin 0 sin(<t> - 4>) - z t ' cos 0. (52) 

The optical path length of the transmitted ray with 
respect to that of the reference ray is then 

Strain = <*o + ns 0l + 01 - a* - pfl = z 0 ' + ns 0l 
- Ar{ cos sin 0 cos(<$ - 4)) - Brf 
x sin £ { ’ sin 0 sin(0 - 4>) - zf cos 0, (53) 

independent of whether is on the shadowed or lit 
side of the spheroid. 

Each ray transmitted through a spherical particle 
participates in two focusing caustics 32-34 : (i) a spheri- 

cal aberration cusp of revolution, which is also known 
as the tangential caustic, and (ii) an axial spike caustic 
on the cusp axis, which is also known as the sagittal 
caustic. The cusp of revolution has its apex at the 
paraxial focal point and it extends through the sphere 
surface to its interior. 34 The axial spike caustic 
meets the cusp of revolution at the paraxial focal 
point. Using van de Hulst's rule for equating the 
number of -tt/2 jumps in the phase of a ray to the 
number of caustic participations 35 and including an- 
other factor of — ff/2 to compensate for the overall 
factor of((0jin Eq. (3), we obtain 

2t r 3-rr 

^ ^trana \ ^trons ~ « (84) 

A. I 

p<V^v\f^W> . rt'i'VV CAM/U. CtW £U/> jvi rtk l-'J u l" 

as the phase of the transmitted ray/ 

For an arbitrary spheroid eccentricity and orienta- 
tion the situation is considerably more complicated. 
Nye's analysis of the one-internai-reflection rays for 
side-on incidence 19 gives a hint as to the complexity of 
the focusing caustics of spheriods. We have not yet 
fully solved the caustic problem for transmission 
through an arbitrarily oriented spheroid. But in the 
remainder of this section the complete solution for 
end-on incidence is given first and then a few com- 
ments concerning the case of arbitrary orientation 
are made. In Fig. 3 we show the evolution of the 
transmission caustics as b/a is increased from 1.0 for 
the specific example of n = 1.333. Figure 3(a) corre- 
sponds to transmission through a sphere, where the 
cusp caustic of revolution points outward from the 
sphere. Figure 3(b) corresponds to b/a = 1.33. 
The cusp caustic begins to retract into the spheroid 
while continuing to point outward. Figure 3(c) corre- 
sponds to 
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— _ al . thirci focusing caustic for this case. For b 2 /a 2 only 

a 2 ~ n 2 - 1 ’ ^5) slightly larger than n/(n - 1), the rays with small r Q ' 

exit the spheroid with a negative scattering angle, 
where all the incident rays focus at the point 36 ’ 37 As r fl ' is increased, the scattering angle further 

decreases, reaches a relative minimum, 38 and then 
z> * (d 2 ~ u 2 ) l/2 , (56) increases, reaching© = 0° when 


independent of r 0 ' and §>'. In Figs. 3(d)— 3(h) the 
interior cusp caustic of revolution points inward 
toward the center of the spheroid. Thus the se- 
quence in Figs. 3(a)-3(d) may be interpreted as 
turning the cusp of revolution inside out through its 
contraction to a point focus in Fig. 3(c). 

Figure 3(e) corresponds to 33 

b 2 n 


where the paraxial rays cross the z axis at the origin 
and exit the spheroid with the scattering angle 0 = 
0°, giving a forward glory. The forward giory is a 


<7o 



b 2 /a 2 a 2 


(58) 


For yet larger values of r 0 \ the scattering angle then 
becomes positive. The rays between r 0 ' = 0 and the 
value of r 0 ' corresponding to the relative minimum of 
0 (i.e., r 0 ,fnm ) form a third and fourth focusingcaustic. 
These rays-are a new cusp of revolution far outside 
the spheroid and pointing inward and a new axial 
spike caustic. The new cusp caustic evolves into a 
far-zone transmission rainbow 33 at the minimum 
scattering angle 0 min corresponding to r 0 ' m,n . For r 0 ' 
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Fig. 3. Rays transmitted through a spheroid for a * 1.333. (a) 6/a* 1 corresponding to a sphere. The caustics are a cusp of revolution 
pointing outward and an axial spike caustic. (b)6/a =• 1.33. The cusp starts to retract into the spheroid. (c)6/a * 1.5123. All the rays 
focus at a singie point, (d) 6/a =3 L.67. The caustics are a cusp of revolution pointing inward and an axial spike caustic, (e) 6/a = 2.0. 
Paraxial rays have a scattering angle of 9 =* O’ corresponding to a forward glory, (f) 6/a * 2.25. The paraxial rays form a butterfly of 
revolution caustic and a second axial spike caustic, (g) Butterfly caustic of (0 magnified by a factor of 11.67. The caustic begins at 
location 1 and then continues in order to the locations 2, 3, 4, 5. (h) 6/ a 13 2.45. The second*cusp caustic points outward. As a function 

ot r 0 ' there are two critical angles for total intemai reflection. 


between r 0 ' min and the value implicitly given in Eq. 
(58) only the new axial spike caustic is formed. It 
extends from beyond the second cusp point out to 
infinity, where it evolves into the forward glory. 

As b 2 / a 2 is further increased, the new inward* 
pointing cusp also begins to turn itself inside out by 
progressing through a butterfly of revolution caus- 
tic, 39 as is shown in Fig. 3(f) and greatly magnified in 
Fig. 3(g). Finally, if b 2 /a 2 is sufficiently large, the 
evolution of the butterfly caustic reaches completion, 
and it becomes an outward-pointing cusp of revolu- 
tion, as in Fig. 3(h). Although Figs. 3(a)-3(h) de- 
scribe n - 1.333. the same evolution of the caustics 
was observed for all the other refractive indices 
examined for the end-on spheroid geometry. 


As a result, depending on the values ofrc, b 2 /a 2 , and 
r Q ' for end-on incidence, a transmitted ray partici- 
pates in either 2, 3, or 4 focusing caustics and 
acquires a transmitted phase of either -3-/2, -2 
or — 5 tt/ 2, respectively [including an additional factor 
of — — /2 to compensate for the overall factor of i in 
Eq. (3)j. The number of caustic participations is 
illustrated in the phase diagram of Fig. 4. The case 
of an oblate end-on spheroid presents no special 
problem because the rays transmitted through it 
resemble the rays transmitted through a thick lens. 
Each ray participates in only 2 caustics, a spherical 
aberration cusp of revolution and its associated axial 
spike caustic. 33 

If the z axis of the spheroid is now tilted with 
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Fig. 4. Phase-space diagram for the caustics produced by the 
transmitted rays when r 0 ' and b/a are varied while the refractive 
index n is held fixed. In region A there are two caustics, a cusp of 
revolution pointing outward and an axiai spike. Along the line ota 
corresponding to Eq. (55) these caustics contract to a point 
focus. In region B there are two caustics, a cusp of revolution 
pointing inward and an axial spike. In region C there arc three 
caustics, a cusp of revolution pointing inward and two axiai 
spikes. In regions D, E, F, there are four caustics. In D they are 
two cusps pointing inward and two axiai spikes. In E they are an 
inward-pointing cusp, a butterfly, and two axial spikes. In F they 
are an inward-pointing and an outward-pointing cusp, and two 
axial spikes. The line 3(J corresponds to the forward glory given 
by Eq. (58), and the line Is the transmission rainbow. The 
cross-hatched regions denote the absence of transmitted rays as a 
result of total internal reflection as in Eq. (33). 


respect to the lab z axis for b 2 /a 2 < n/(n - 1), the 
rotational symmetry that produced the degeneracy of 
the axial spike caustic is lost and it deforms into a 
four-cusped astroid caustic. 40 The cusp caustic also 
loses its rotational symmetry. 41 The combination of 
the distorted cusp and astroid has an astigmatic 
focusing character, and sections through the compos- 
ite caustic resemble the evolution shown in appendix 
2 of Ref. 39. For b 2 /a 2 > n/{n - 1) the new cusp 
and axial spike caustics deform into a second compos- 
ite astigmatic focusing caustic as in appendix 2 of Ref. 
39. For yet larger values of b/a, the butterfly of 
revolution caustic evolves into a complicated struc- 
ture that we do not yet fully understand. But in any 
event, as long as the spheroid eccentricity satisfies 
b 2 /a 2 < ri/(n - 1), each ray participates in two 
focusing caustics, and the transmitted phase shift for 
an arbitrary spheroid orientation is given by Eq. (54). 

C. Polarization of the Transmitted Electric Field 

The incident plane wave is polarized so that its 
electric field makes an angle x %v ith the x axis. The 
polarization vector of the incident plane wave with 
respect to the x’y'z' rotated lab axes is then (equation 
48 of part I) 

l xnc = cos(x - ®)u z ' -f- sin(x - 4>)d/. (59) 

In part I, the unit vectors in the TE and TM 


TE tl 



Fig. 5. TE and TM polarization directions for the reflected ray 
and the transmitted ray of Eq. (60) and Eqs. (61H64), respectively. 


polarization directions were taken to be (equations 49 
and 50 of part I) 


TE; 


nc 


sin 0j nc 


— A.„c x{m X k itte ) 
1 E mP = * 


sin 0 in 


(60) 


When considering the transmitted ray we define the 
TE and TM polarization unit vectors to be the 
negative of Eq. (60). The reason for this is as 
follows. For scattering by a sphere, an observer in 
the 0, direction will intercept a reflected and a 
transmitted ray that originated on opposite sides of ^ ^ 
the sphere. This is shown in Fig. 5. These re- ^ 
fleeted and transmitted fields are added together to 
form the total field at the observer. The addition is 
most simply performed if each of the two individual 
fields is decomposed into components in the same 
directions; i.e., the outgoing TE and TM directions for 
reflection are the same as the outgoing TE and TM 
directions for transmission. This is ensured for ray 
scattering by a sphere if the TE and TM directions for 
reflection are given by Eq. (60) and the TE and TM 
directions for transmission are given by 


TE 


t0 


'fetO X rn£ 
sin 9j 0 


— . k tQ X (k fQ X rh Q ) 

™:0 = — a 

sin 0 tO 


(61) 


TT? 

1 


k t Q x 
sin 9; 0 
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trans 



tm !0 = 


kto x (^fo x rho) 
sin 9f 0 


— k n x fi*' 

TE “ = ImC ’ 

- — k a x (k a x fi,' 

TM U = 1 1 1 


sin 9 :1 


— - k n x n/ 

TE ;1 = 


TM (l = 


sin 0 n 

x (*n * ^i') 
sin 9 f , 


(62) 


(63) 


(64) 


as is shown in Fig. 5. 

For scattering by a spheroid, the rays reflected and 
transmitted in the 0, direction again in general 
originate on opposite ends of the lit side of the 
spheroid. As a result,' we use the TE and TM 
polarization directions of Eqs. (61H64) for spheroid 
scattering as well. For an arbitrary spheroid orienta- 
tion with t) Q t] L , the two TE unit vectors TE f0 and 
TE tb which describe the ray inside the spheroid 
referenced with respect to the 0 interface and the 1 
interface, respectively, do not coincide. This is the 
source of the first of the cross-polarization effects in 
the transmitted intensity. We derive this effect as 
follows. The incident-ray polarization vector of Eq. 
(59) may be decomposed into the TE and TM compo- 
nents of Eq. (61), giving 


- £inc * cos 7 oTE ! 0 + sin 7 0 TM e0 , (65) 


where 


7o - X - 4> - ^0 - 9 • (66) 

After transmission into the spheroid, the polarization 
vector of the ray becomes 


£f 0 ~~ ^te' 3 cos 7 qTE £ q -r sin YoTM fQ (67) 

where t TE °{Q iQ ) and are the Fresnel coeffi- 

cients for transmission corresponding to the angle of 

incidencej}^. As mentioned above, the TE, 0l TM (0 

and the TE tl , TM n polarization vectors are rotated 
with respect to each other because the 0-interface and 
the 1-interface planes of incidence do not coincide. 
Specifically, by substituting Eqs. (11), (15), and (27) 
into Eqs. (62) and (63) we obtain 

TE f0 = cos A 01 TE a 4* sin d 0l TM tlJ 

TM f0 — -sin Aq^TE^i t cos AqjTM, (68) 

where 

v si n^ t sin(T] 0 -p : ) 

^an ig, , , i- 1— . . - 

sm n' l0 cos 4* cos M' l0 sin 4 / 1 cos(ti 0 - 


E 



Fig. 6. Contributions to the transmitted electric fieid. The 
contributions proportional to cos A oi are the polarization terms 
that describe transmission by a sphere or end-on spheroid. The 
contributions proportional to sin a 0I are the first type of cross- 
polarization terms. 


if the ray exits the spheroid on its shadowed side. 

If it exits on the lit side, the factor of cos I', in Eq. (69) 
is replaced by -cos 'P l . The polarization vector of 
the ray incident at the point of exit on either the 
shadowed or lit side of the spheroid is then 

«;i = ('te 0 cos 7 o cos A 0I - C TM ° sin y a sin d 01 )fE a 

+ (t TE ° cos 7 0 sin A 01 + i TM \ sin y 0 cos A 01 )TM ;I . 0 

(70) 

Finally, the polarization vector of the outgoing ray 
after it exits the spheroid is , 

^ . J, A wtvua 

<=!! =(t TE l, t TE I cos 7 0 cosii 0I - /f TM °i T ^sin 7 0 sin^ 0 ,)fE fl 




Ute^tm 1 cos 7 0 sin d 01 

7o cos A 01 )TM,,, 


(71) 

where iTE^(Oii) and are the Fresnel coeffi- 

dents for transmission corresponding to the angle of 
incidence 9^. Again, 9 tl must be less than the critical 
angle for total internal reflection in order for the ray 
to exit the spheroid. 

4. Cross-Polarization Effects r ' 

A pictorial representation of Eq. (71) /b e lW is given 
in Fig. 6. For scattering by a sphere, the only 
possibilities are TE polarization in producing TE 
polarization out and TM polarization in producing 
TM polarization out. For scattering by an arbi- 
trarily oriented spheroid, these two possibilities are 
weighted by cos A 0l . The cross-polarization possibili- 
ties of TE polarization in producing TM polarization 






r 
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flected-ray TE and TM polarization basis while diffrac- 
tion is diagonal in any polarization basis. For trans- 
mission, the rightmost matrix converts the incident 
reflected-ray-polarization basis to the incident trans- 
mitted-ray polarization basis, and the leftmost matrix 
converts the final transmitted-ray polarization basis 
back to the final reflected-ray polarization basis. 
Transmission at each interface is diagonal in the 
transmitted-ray TE and TM polarization basis, and 
the middle matrix converts the O-mterface transmit- 
ted-ray polarization basis to the 1-interface transmit- 
ted-ray polarization basis. Explicitly, the compo- 
nents of the scattering matrix are 

Si(0, <t>) = S diff + exp(i5 rt f)r TE + S trans 

x exp(: 5 trans )(ixs 0 f TE l cos cos Aq! cos fl 
— £xe^tm^ sin sin Aq^ cos I"! ^tm ^te 
x cos A R1 sin A 0l sin £1 - si n ^si 

x cos A 0 . sin fl), {?*?) 

S->(Q, <J>) = S dif r + S„( exp(£5 r rf)r TM + ^trana 


Re = TE C0S 

R\f = *S ref r TM sin y, 

Trr = S trans ^ T £°f TE L COS y 0 , 

T £M = S trans ^E 0f TM l C0S 7o> 

T\tE = 5 tr3ns ^TM°^TE I 7o» y 

Tmm = S trans /xM 0iJ TM l sin 70* (8T) 


we obtain 


ref+trans(®> 


gp 2 1 

2^0 k-R 2 


[Sdifr 2 + Re- + R« 2 


r 



+ ('R EE~ + Tvf«')cOS* A 01 

+ 2(i?rZ?£ + RMcos 5 re f 
+ 2 (T ee D e + T MM D M )cos A SI cos A 0l cos ^ trans 
+ 2 (T EE f?r 4 “ TvfA/i ?.^) 003 
x cos A 0l cos(5 trena - Sfrf) 


x exp(io tnm9 )(^ TM 0 fTM l COS A ffl cos A 01 cos fl 


+ dW + Taft sin 2 A 


oi 


— ^tm^te 1 sin A rtl sin A 0l cos fl - *xe°*tm 1 

x cos A H1 sin A 0l sin fl - ^te%e 1 s ^ n A R\ 

x cos A 0l sin fl), (38) 

SO 

^3(®» “ Stran* exp(i5 tran3 )(i T M°^TM I cos &ri cos ^oi 

X sin fi 4* i^TE%M l cos ^K 1 sin -Iqi cos ^ 

4" ^XE^TE^ sin l\ R i COS Aqj cos 0 

— ^XM°^TE l sin ^Rl s ^ n -Iqi Sin ^)» C**®) 

o 1 

S 4 (0, <&) = S tnuis exp(i5*„i ng )(— ixs°fxE l cos cos ^oi 
x. sin fl — ^xm^te^ cos A^j sin Aq^ cos fl 

— ^TM^TM 1 s ^ n COS A 01 COS fl 

4- ^xe^tm 1 sin A/ji sin A 0l sin fl). (Stf) 

For a spheroid that is tilted away from end-on 
incidence by only a small amount, the first transmit- 
ted terms of S\ and S<> are comparable to the dif- 
fracted and the reflected terms, and the last three 
transmitted terms are second-order corrections; i.e., 
they have two terms in the sine of a small angle. 
Similarly the first three terms of S 3 and S 4 are first 
order in the sine of a small angle, and the last term is 
of the third order. 

Finally we obtain the ray-theory scattered intensity 
by multiplying Eq. (72)15y its complex conjugate after 
inserting Eq. (.76]^ The resulting expression is rather 
long. But using the shorthand notation, 

De = S dirr cosy, 

D\t - Sdiff sin y, 


4- 2{T mm T km - T EE 7\, E )cos A 0l sin A 0l 
4- 2(T EE i?.u - T MM i?r)sin A rtl 
x cos A 01 cos(5 trans - Sprf) 

4“ 2{Tz:mR\i “ T M R s )ca* 
x sin A 01 cos(5 trans - b n( ) 

~ 2(T iV f E i?v + T £M R E )sin A 
x sin A 01 cos(5 trnnj - 5„f) 

4“ 2 i^T ee^ st T MM D s)siii A/ji cos A 0l cos 5. rans 

4- 2{T em D m — T^gD?) cos Aft i sin A 0l cos 5> rans 
— 2(T r J vf E i).\f 4- r^DpJsin A#i sin A 0l cos ^transj- L. 



The various terms in Eq. (82) haye th^Jbllowing 
physical interpretations. The second line j^the indi- 
vidual diffracted, reflected, and transmitted intensi- 
ties for scattering by a sphere or an end-on spheroid. 
The reflected and transmitted intensities in the 0, 
direction depend on the polarization sta te [E or MXsL 
the incident beam. The th4rd~~an d fourth - Imesare 
the diffraction-reflection interference, the diffraction- 
transmission interference, and the reflectiorb^ trans.- 
mission interference for scattering by a sphere or 
end-on spheroid. The fifth lin e al s the cross-polariza- 
tion contribution to the transmitted intensity, whj^h 
is due solely to the rotation of the TE f0 and TE tl 
polarization directions. For scattering in the near- 
forward direction by a spheroid tilted a small amount 
from end-on incidence these two terms are second- 
order corrections because sin A 01 is small and because 
the TE and TM Fresnel coefficients are nearly equal 


-re--*' 

“ 4vil 
•j ( v-l 'S\ ’ 

j i 
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for small angles of incidence. 44 The -sixth through 
- — lines are the cross-polarization contribution to 
the reflection-transmissicfo interference. The mn4h 
through _^ftven4h lines are the cross-polarization 
contribution to the diffraction-transmission interfer- 
ence. For scattering in the near-forward direction 
by a spheroid tilted a small amount from end-on 
incidence these are also second-order corrections, 
either because of sine-squared factors or because of a 
single sine factor and the near equality of the TE and 
TM Fresnel coefficients. Because diffraction and 
reflection are both diagonal in the reflected-ray basis 
(see Eq. (76jfthere is no cross-polarized diffraction- 
reflection interference. The cross-polarized contribu- 
tions to Eq. (82J uecome comparable to the second 
^__thrnug^ lines for a highly eccentric spheroid 

" tilted substantially away from end-on incidence and 
for rays with large angles of incidence at the points of 
entrance to and exit from the spheroid. 

5. Computation of the Scattered Intensity 

It has already been shown in part I that the diffracted 
plus reflected electric field is directly expressible in 
terms of 0 and But because of the complexity of 
both the spheroid shape and the refraction geometry, 
the transmitted electric field is not directly express- 
ible in terms of 0 and <t>. This suggests the following 
strategy for computing the ray-theory intensity ofEq. 

We first generate a dense grid of r 0 ' and £q' 
values, i.e., Ar 0 ' = 0.005 and = 1°. Then for 
each r 0 \ 5 0 ', we compute the scattering angles 0 and 
<t>. We then numerically perform the derivatives in 
Eq. (46) and obtain the magnitude, phase, and polar- 
ization direction of the transmitted ray. For each 
r 0 ', §o' also check whether the ray is totally 
internally reflected by the use ot Eq. (33). If it is not 
totally internally reflected, we test whether it exits on 
the shadowed or the lit side of the spheroid with Eqs. 
(25) and (26). Knowing © and for each transmit- 
ted ray, we then calculate the reflected and diffracted 
electric fields for those scattering angles, and finally 
we compute the scattered intensity of Eq. (&2JAH 
In order to test our numerical procedure, first b = a 
was set in the computer program, and the results 
were compared against ray scattering by a sphere as 
given by the analytical formulas in Refs. 3 and 4. 
The results matched exactly. For scattering by a 
sphere or an end-on spheroid, Eq. (46) reduces to 


The agreement between our computer program and 
the numerical implementation of the analytical formu- 
las for sphere scattering in Refs. 3 and 4 verifies that 
we have chosen our grid of r 0 ' values to be fine enough^ 
to perform the numerical derivative in Eq. (83r 
accurately. Our computed result for b — a — 10.07 
pirn also closely approximates^Mi^ theory, as is shown 
in Fig. 3. Also shown in Fig. 3 is the generalized 
eikonal approximation. 45 The eikonai approxima- 
tion has been shown to be an accurate approximation 


in the short-wavelength limit. It has the feature 
that it slightly underestimates the peak intensity of 
the reflection-transmission interference structure for 
20° ^ 0 S 50° forn = 1.33. 

We next considered the end-on spheroid case \ = 
0.6328 p.m, n = 1.33, b/a = 1.5, a - 10. 0 am, x = 90°, 
and 9 = 0° and compared our results with those in 
figure 8 of Ref. 21, in which the ray-theory intensity 
was calculated with another method. Again the 
results matched exactly. 46 We also compared our 
results with the generalized eikonal approximation 
for a = 10.07 pun. 47 The comparison is shown in Fig. { 
9(a). For 0 < 30° ray theory and the generalized 
eikonal approximation give similar results, with the 
generalized eikonal approximation again slightly un- 
derestimating the reflection-transmission interfer- 
ence for 0 > 10°. The critical angle for total inter- 
nal reflection occurs at a scattering angle of 
approximately 0 = 38°. Thus ray theory is expected 
to be inaccurate for 0 > 33° because of our neglect of 
the Fock transition 48 in the transmitted electric field. 
This neglect may be the cause of ray theory becoming 
out of phase with the eikonal approximation for 0 > 
25°. 

In Fig. 9(b), ray theory is compared with the 
generalized eikonal approximation for side-on inci- 
dence with \ = 0.6328 jum, n = 1.333, b/a - 1.5, a = 
10.071 p.m, 9 = 90°, <$> = 0°, and x = 90° for scattering 
in the = 0° plane. The cross section of the 
spheroid in the xz plane (i.e., — 0°) resembles an 

oblate spheroid with b/a = 0.667 in the end-on . 
incidence configuration. For this configuration and 
c f> = 0°, each transmitted ray participates in two 
focusing caustics. Again ray theory and the general- 
ized eikonal approximation give similar results, with 
the generalized eikonal approximation slightly under- 
estimating the reflection-transmission interference 
for 15° < 0 < 40°. The critical angle for total 
internal reflection occurs at a scattering angle of 
approximately 0 = 55°. Thus our neglect of the 
Fock transition causes ray theory to be inaccurate for 
0 > 50°. 

In Fig. 10 ray theory is compared with the exact i 
solution of the plane-wave-spheroid problem with 
the method of Ref. 49. The comparison was made 
for end-on [Fig. 10(a)) and side-on [Fig. 10(b)) inci- 
dence for X = 0.6328 pLm, n = 1.333, b/a — 1.5, a — 
3.021 p.m, and x “ 90° for scattering in the T = 0° 
plane. This corresponds to a spheroid-size param- 
eter of 2— /\ = 30.0, which is at the lower end of the 
region of applicability of ray theory. In Fig. 10(a) the 
comparison is good for 0 < 30°, and in Fig. 10(b) 
the comparison is good for 0 < 50°. We consider the 
general agreement between our results and the gener- 
aiized eikonal model for 2 rra/\ = 100 and the exact 
W solution for 27 ra/\ - 30 as an additional check of the 
correctness of our method. 

It is of great interest to compute the cross- 
polarization contributions to the scattered intensity 
of Eq. Unfortunately, this cannot be done 

reliably until a complete solution to the problem of 
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Fig. 3. Intensity as a function of the scattering angle 0 with <P * 0* for a plane wave with \ = 0.632S p. m and x 3 90“ incident upon a 
sphere with a =» 10.071 pirn and n - 1.333. The solid curve is th<M!3theory result, the dashed curve is Eq. (827, and the dotted curve is the 
generalized eikonal approximation of Ref. 45. (. 
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the number of caustic participations of an arbitrary 
transmitted ray is obtained. This problem is cur- 
rently under study. It is also of great interest to 
compare our results with experimental data. 
Although' much experimental spheroid scattering data 
concerning the generalized rainbow in the backward 
hemisphere exists, 16,50 “ 53 experimental spheroid scat- 
tering data in the near-forward direction appear to be 
sparse. 54 

Last, it has been pointed out by a number of 
authors 3 * 11 * 12 that both the T-matrix solution and the 
spheroidal wave-function solution of the problem of 
scattering a plane wave by an arbitrarily oriented 
spheroid possess numerical difficulties for prolate 
spheroids with large b/a ratios.. In ray theory, we 
have found that a number of alternatives occur for 
rays transmitted through large b/a prolate spheroids. 
These include the point of exit of the ray being on 
either the shadowed or the lit side of the spheroid, a 
complicated caustic structure that dictates the phase 
of the transmitted electric field, and the existence of 
more than one Fock transition of the transmitted 
electric field as a function of scattering angle [see Fig. 
3(h), for example]. It would be of interest to deter- 
mine whether any of these complexities for eccentric 
prolate spheroids are related to the difficulties in 
computation for the associated wave-scattering prob- 
lem. 

The work was supported in part by the National 
Aeronautics and Space Administration gTant NCC 
3-204. I wish to thank T. W. Chen and L. Yang of 
New Mexico State University for permitting me to 


quote some of their results before publication. I also 
thank J. P. Barton of the University of Nebraska, 
Lincoln, for providing me with the exact scattering 
data quoted in Fig. 10. 
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